Normal metal - superconductor tunnel junction as a Brownian refrigerator 
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Thermal noise generated by a hot resistor (resistance R) can, under proper conditions, catalyze 
heat removal from a cold normal metal (N) in contact with a superconductor (S) via a tunnel 
barrier. Such a NIS junction acts as Maxwell's demon, rectifying the heat flow. Upon reversal of 
the temperature gradient between the resistor and the junction the heat fluxes are reversed: this 
presents a regime which is not accessible in an ordinary voltage-biased NIS structure. We obtain 
analytical results for the cooling performance in an idealized high impedance environment, and 
perform numerical calculations for general R. We conclude by assessing the experimental feasibility 
of the proposed effect. 
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In 1867, Maxwell suggested a demon that attempts 
to violate the second law of thermodynamics [lj. The 
demon acts between two containers A and B, initially at 
the same temperature; it exclusively allows hot particles 
to pass from container A to container B, and cold ones 
from B to A. This process would lead to a decrease of 
entropy if the system were isolated, and could then be 
used for useful work. Ever since, this thought experiment 
has intrigued physicists, see, e.g., Refs. [10, [111. Yet 
the demon needs to exchange energy with the containers 
in order to function properly. Thereby the entropy of the 
whole system (including the demon) is always increasing, 
rendering thermodynamics intact. 

In this Letter we present a particularly illustrative ex- 
ample of Maxwell's demon, which can be realized ex- 
perimentally in a straightforward way. Our system is a 
Brownian refrigerator [5| in close analogy to Brownian 
motors and thermal ratchets [1, 0, H, 0, G3- It con- 
veys heat unidirectionally in response to random noise. 
Specifically, we consider a tunnel junction between a nor- 
mal metal and a superconductor (NIS junction) subject 
to the thermal noise of a resistor at temperature Tr, 
see Fig. [TJ The temperatures of the electrodes N and 
S are Tn and Ts, respectively. The capacitance C con- 
sists of that of the junction itself and the surrounding 
circuit. The resistor and the junction can be connected 
by superconducting lines which efficiently suppress elec- 
tronic thermal conductance, and thus enable us to dis- 
cuss solel y t he photonic heat exchange via the lines as 
in Refs. [IJ, [l2|] . The N side can be connected to the 
superconducting line via a metal-to-metal SN contact, 
which provides perfect electrical transmission but, due 
to Andreev reflection, prevents heat flow [H, fl4| . We 
note already here that the presented NIS structure can 
be replaced by a more conventional symmetric SINIS de- 
vice with two tunnel junctions back-to-back. This will be 
discussed towards the end of the present Letter. We also 
note that the idealized circuit of Fig. [TJ (a) as well as the 
thermal model presented schematically in Fig. [TJ (b) and 
mathematically below are not unrealistic for a practical 
on-chip device. 

The heat balance between the resistor and the NIS 
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FIG. 1: Schematic presentation of the system. In (a) we show 
the electrical diagram of the resistor at temperature Tr, and 
the tunnel junction. The parallel capacitance C includes that 
of the junction and a possible shunt capacitor. In (b) we show 
the thermal diagram of the system. The NIS tunnel junction 
acts as a Brownian heat engine, or as Maxwell's demon (MD) 
between N and S. Qn and Qs are the heat fluxes out from 
N and S, respectively, and P e xt denotes the external power 
needed to create the temperature bias of the resistor. The 
thermal resistances indicate phonon coupling of the electrical 
subsystems. 



structure can be described on the level of a single elec- 
tron tunnelling event between N and S, accompanied by 
an exchange of energy with the junction's resistive envi- 
ronment. Formally this can be done using the so-called 
P(E) theory developed for a tunnel junction embedded 
in an electromagnetic environment [TH fl6| . Assuming 
that for both electrodes i = 1,2 the (normalized) den- 
sity of states rii (E) is symmetric around the Fermi level 
E = and that the corresponding energy distributions 
satisfy fi(E) = 1 — fi(—E), we find that the heat flux 
out of 1 upon tunnelling from electrode 1 to 2 is given by 
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Here, P(E' — E) is the probability density of emitting 
energy E' — E to the environment in a tunnelling process 
from 1 to 2. It can be calculated for a particular envi- 
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ronment at a temperature Tr, once the dissipative part 
5RcZ(cj) of its impedance at frequency uj/2ir is known 
[lf|. The theory is perturbative in tunnel conductance, 
therefore the tunnel resistance i?x should be of the or- 
der of the resistance quantum, i?x = h/e 2 , or higher. 
Here, we have neglected the Joule dissipation in the nor- 
mal metal itself, which is usually justified because N can 
have a very low resistance as compared to that of the 
tunnel junction. 

As a warm-up exercise let us look at the simplest 
system, where the tunnel junction is of type NIN, i.e., 
both sides are normal metals. Then rii(E) = 1 for both 
electrodes; we furthermore assume them to be kept at 
equal temperature Tjv and characterized by equilibrium 
Fermi distributions /, (E) = (1 + e E/k B T N yi_ Thig gx _ 
ample closely resembles the original set-up of Johnson 
and Nyquist [l7l |. where thermal noise acts between two 
resistors. At high temperatures, ^bTrn 3> h(RC)^ 1 , we 
then obtain, using Eq. {TJ, the total heat flux out of the 
junction as Q = Q x + Q 2 ^ k B (T N - T R )(i? T C)- 1 . This 
power is shared equally by both electrodes due to symme- 
try. The result for Q is consistent with the expectation 
that heat flows from hot to cold. 

The symmetry of the previous example is vitally bro- 
ken in the NIS junction that we focus on here. In par- 
ticular, we have n\{E) — 1 for N as before, but the BCS 
density of states in S, n 2 (E) — for \E\ < A and n 2 (E) = 
\E\j\jE 2 — A 2 for \E\ > A, makes this system behave as 
an energy selective entity in the spirit of Maxwell's de- 
mon. Here, A is the energy gap of the superconductor. 
We again assume standard Fermi distributions, but pos- 
sibly at different temperatures: fi(E) = (1 + e - E/ ' CBTN )~ 1 
for N and f 2 (E) = (1 + e E / k ^ Ts )- x for S. 

Let us first discuss a particular limit where illustra- 
tive analytical results can be obtained. We consider 
a very resistive environment, such that 7r^- ^ h 3> 

2 

1, where Eq = §q is the charging energy. Then 
P(E) assumes a simple Gaussian form: P(E) = 

(27rg)~ 1 / 2 exp(— ^ E ~2^ )i where the width is given by 
a = 2fc B TR£'c [HI]. We now assume that this width is 
small compared to the superconducting gap A. The in- 
tegration over energy in {1} then involves E, E' > A ^> 
&b7n,s and we may approximate the Fermi distributions 
by their Boltzmann-like exponential tails. In this limit, 
setting further Tg = Tn for simplicity, we have 
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One can similarly estimate the heat flux from the super- 
conductor: 
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The expressions ^ and vanish when Tr = T^, as 
expected. The non-trivial result is that on heating the 



environment to temperatures Tr > Tn, N tends to cool 
down, i.e., Qn > 0, and S is heated, Qs < 0. On the 
other hand, for Tr < Tn the heat flux is reversed: S 
tends to cool down and N to warm up. Also in this 
regime the heat flux between the junction and the resis- 
tor is unevenly distributed among N and S. Such cooling 
of S never occurs in a conventional voltage biased NIS 
refrigerator [l4| . 

If one employs ([2]) to find the optimum Tr/Tn where 
the cooling power of N is maximal, one finds Tr/Tn — 
A/2E C for E c /A < 1. Although this is the right or- 
der of magnitude, the numerical results presented below 
show that the actual value of the ratio Tr/Tn is approx- 
imately twice higher. The approximation used above for 
the Fermi functions leads to a counterbalance of the cool- 
ing effect as the Boltzmann tails grow exponentially at 
negative energies. A better estimate can be obtained by 
retaining the full Fermi function fi(E') and linearizing 
the exponent of the Gaussian P(E' — E) around E' = 
where /i steeply drops. We then arrive at 
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This expression predicts that for &bTn, Eq <C A the op- 
timal point of cooling indeed lies at Tr/Tn — A/Eq. In 
its range of validity this optimum as well as the overall 
behavior described by Eq. (j4|) are consistent with the 
numerical results, as will be seen below. 

In order to picture the characteristics of the system 
in a dissipative environment of arbitrary resistance R, 
we have performed numerical calculations. The usual 
environment theory of single-electron tunneling was em- 
ployed here also. In all the numerical calculations the 
ideal junction was coupled to a parallel RC circuit. In 
Fig. [2] we show a comparison of the exact numerical re- 
sults for R i?K and the analytical approximations of 
Eqs. © and We see in the main frame that at low 
temperatures, in this case at /cbTn/A = 0.03, the approx- 
imation (j4|) works quite well over a broad range of tem- 
perature biases, whereas the Boltzmann approximation 
((2]) fails at high values of Tr/Tn- Yet, around zero tem- 
perature bias the latter approximation works perfectly as 
demonstrated by the inset of Fig. [2l 

A set of numerical results under representative condi- 
tions is collected in Fig. [3l In (a) we see the influence of 
Ec on the performance of the system. Here the Gaussian 
approximation of P(E) was used. The maximum cool- 
ing occurs indeed at Tr/Tn — A/Eq, and the value at 
the maximum grows a little with increasing A/Eq. In 
(b) the Gaussian approximation was abandoned. Under 
reduced R the P(E) function transforms from a broad 
Gaussian (width 2/cbTr-Ec) centered around Ec towards 
a delta- function around E = 0. According to Eq. {J) this 
evolution weakens the refrigeration effect. In Fig. [3] (b) 
we see that it is indeed essential to have a relatively large 
R in order to sustain the effect, although some cooling 
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FIG. 2: Calculated cooling power Qn at ksT^/A = 0.03 as 
a function of the temperature bias Tr/Tn. Here we have 
assumed that R 2> Rk, Ts = Tn, and A/Ec = 5. In the 
main frame the solid line is the result of the exact numerical 
calculation, and the dashed one is the analytic expression of 
Eq. @ based on linearizing the exponent of the Gaussian 
distribution over the relevant energy interval. The dotted line 
diving to negative values around Tr/Tn — 3 is the result of 
Eq. ([2]). This approximation works well, contrary to Eq. Q, 
at small temperature biases, as shown by the inset: the solid 
line is the exact result and the approximation @ is shown by 
open dots. 



can be observed also down to R/Rk ~ 0.1. We also see 
that with high environment resistances the cooling char- 
acteristics approach the result of a Gaussian environment 
as should be the case. In (c) and (d) we show the quan- 
titative results for the heat flux Qn and the optimum 
temperature bias, respectively, as a function of Tjy in a 
highly resistive environment. 

Although the discussion above demonstrates counter- 
intuitive heat fluxes in the system, it is straightforward 
to verify that they do not contradict the second law of 
thermodynamics. For instance, consider the resistor and 
the NIS junction to form an isolated system such that, 
referring to Fig. Q] (b), P cxt = and that couplings to 
the phonon bath vanish. Initially brought to tempera- 
tures Tr and Tn (= Tg for simplicity), the system then 
starts to relax towards a common temperature. For a 
small temperature difference AT = Tr — Tn, we have 
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cording to Eqs. ((2|) and ([3]). This is the heat flux between 
the junction and the resistor, and it is always directed 
from hot to cold thereby ensuring positive entropy pro- 
duction. 

Some of the results above can be obtained approxi- 
mately by straightforward classical estimates. In par- 
ticular, we may consider an ordinary NIS junction 
with large resistance Rt and evaluate its cooling power 
when biased by a fluctuating voltage with vanishing 
mean value and Gaussian variance ((el^) 2 ) = a = 
2&b7r£'c- hi the limit of low frequencies (small 
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FIG. 3: Heat flux out from the normal island. In (a) and 
(b) we assume that £;bTn/A = 0.1 and that the supercon- 
ductor temperature Ts = Tn- The cooling power is plotted 
as a function of the ratio Tr/Tn. In (a) the environment 
resistance is R S> -Rk- Different curves, with maxima from 
left to right, correspond to A/Ec = 1, 2, 5, 10 and 30. In (b) 
A/Ec = 5, and the resistance varies from top to bottom as 
R/Rk = oo,10,2,l,0.5,0.25, and 0.125. In (c) we plot the 
maximum cooling power and in (d) the relative temperature 
of the resistor at this optimum point as functions of the N 
island temperature. In (c) and (d) the curves from top to 
bottom are for A/E c = 30, 10, 5, 2, and 1. 



Eq), we may make a quasi-stationary averaging over 
the fluctuations, such that the expected cooling power 
of N is (Q N ) ^ J p(eV)Q N d(eV), where p(eV) = 
{2n(j)- l/2 e- { - eV)2/{2a) is the distribution of fluctuations 
and Qys(eV) is the cooling power of the NIS junction 
at a static bias voltage V. For \eV\ ^> &bTn, one ob- 
tains Q N (eV) ~ V ^ t Tn/2 (A - | e T/|) e -(A-| e y|)/fcBT N 
[l8|. Performing the averaging yields then (Qn) — 

V2*Afe B T N A -A/fc B T N(1 _ 2 T^E^ )e ^E c /k B T r ^ rc _ 
e z itT v in A ' 



e 2 R T ^ T N A 

suit resembles very closely that of Eq. $2$ , in particular 
for Tr > Tn- Yet the classical result above neglects the 
backflow of heat from the junction to the resistor, which 
is an important contribution especially when Tj> < Tn- 

Next we discuss the possible experimental demonstra- 
tion of the presented effect. The basic concept can be 
implemented for instance in a standard on-chip configu- 
ration by employing common N and S metals like cop- 
per and aluminium, respectively. The resistance can be 
formed using a strip of resistive metal such as chromium 
or a metallic alloy. Figure 2] demonstrates the calculated 
temperature reduction, Tn/Tq, of the N island assum- 
ing that the superconductor is well thermalized at the 
phonon bath temperature Tq. The results are shown 
specifically for the case of aluminium as a superconduc- 
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FIG. 4: Temperature reduction in N in a realistic device at 
the optimized resistor temperature. The resistance is assumed 
to be high, R 3> Rk- The calculation has been performed 
for aluminium as a superconductor with A = 200 jj,eV. The 
other parameters are: A/Ec = 5, E = 1 ■ 10 9 WK~ 5 m~ 3 , and 
Q = 1 • 10 -21 m 3 . From top to bottom the curves correspond 
to R T = 100, 10 and 1 kfi. 



tor (A = 200 jtteV, transition temperature Tc — 1 K). 
We have chosen the ratio A/Eq — 5, which, on one 
hand corresponds to realistic junction parameters, and 
on the other hand allows maximum cooling to occur at 
not excessively high resistor temperatures: the largest 
temperature reduction of almost 40% occurs at T ~ 0.2 
K corresponding to Tr ~ 0.7 K. The specific material pa- 
rameter for the electron-phonon coupling in the normal 
metal was chosen to be £ = 1 • 10 9 WK" 5 m - 3 [3. The 
N island volume of f2 = 1 ■ 10~ 21 m 3 can be achieved by 
a standard process. 

In the analysis above the charging energy Eq is due 
to the parallel connection of the junction capacitance 
and an optional shunt capacitance. We have seen that 
the optimum cooling results do not depend particularly 
strongly on Eq, as long as Ec < A holds. Some im- 
provement in performance can, however, be obtained by 
decreasing Ec, see Fig. [3l because the parallel capaci- 
tance filters out the harmful high frequency tail of the 



noise spectrum. Yet the ratio Tr/Tn ~ A/Eq reflects 
the fact that with large shunting capacitance one needs 
to have a hotter noise source to induce sufficiently strong 
fluctuations. From the practical point of view this is not 
advantageous, at least in a conventional on-chip solution, 
because the high temperature of the environment leads 
to parasitic heating of N via the phonons of the substrate. 
Therefore, in a practical realization the choice of C is a 
trade-off. 

We believe that the presented parallel RC environment 
can be realized almost exactly as long as C is kept small, 
say on the level of few fF arising from the junction itself. 
Then the series inductance of a circuit of sub- 100 fjjn 
dimensions can be neglected, because the corresponding 
LC-frequency remains high as compared to the frequency 
band of thermal radiation at sub-K temperatures. 

Finally, the effect can also be realized in a symmetric 
SINIS configuration, instead of using a single NIS junc- 
tion with an NS Andreev mirror. There the same theo- 
retical analysis can be carried over by replacing the en- 
vironment resistance by an effective resistance i?/4, and 
the capacitance by 2C, where R is the true environment 
resistance and C is the capacitance of one junction. To 
take into account the possible charging effects on the is- 
land, a more involved analysis needs to be performed, 
however. 

In summary, we have discussed a Brownian refrigera- 
tor of electrons, which offers an illustrative example of 
Maxwell's demon in the form of a tunnel junction with 
a superconducting energy gap. Its operation is based on 
building blocks whose characteristics and implementa- 
tion are well known. We expect that it yields a substan- 
tial temperature reduction in a straightforward on-chip 
realization. 
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